IDENTIFICATION OF ANTIBIOTIC TARGETS AND CRITICAL POINTS IN 
METABOLIC NETWORKS BASED ON PATHWAY ANALYSIS 

BACKGROUND OF THE INVENTION 

1 - Field of the Invention 

The present invention generally concerns the identification of 
(i) pathways and (ii) critical points in, and (iii) the generation 
5 of mathematical models of, existing and proposed cellular metabolic 
networks comprised of biochemical reactions or mechanisms with 
genetic or non-genetic associations. 

More specifically, the present invention relates to 
computational methods and systems for the analysis and modeling of 
K> cellular metabolic networks so that, inter alia, potential targets 
L £j in support of the directed development of therapeutic agents and 
^ engineered microbial strains may be identified. 

2 . Description of the Prior Art 
M. 2 . 1 General Background 

15 Within a cell of any organism there are complicated networks 

gg of interacting proteins and enzymes that perform certain chemical 
conversions and transformations. These conversions and 

q transformations -- life processes ultimately lead to the 

M= production of the (i) necessary building blocks (biomass 
20 constituents such as amino acids, nucleotides, phospholipids, etc.) 

and (ii) energy requirements of the cell. Environmental substances 
are processed to meet the demands of a living cell through this the 
cell's network of biochemical reactions. 

These biochemical reaction networks primarily involve the use 
25 of enzymes derived from particular genes whose chromosomal location 
and function have been characterized, as well as enzymes inferred 
to be present based on similarity of their genomic sequence to the 
genomic sequences of enzyme-coding genes in other organisms. There 
is presently, circa 2000, much focus on attempting to model and to 
30 'reconstruct' these networks of a living organism based, primarily, 
on the use of genome sequence information of the organism. 




Meanwhile, the arsenal of reactions that a cell has at its 
disposal dictate the production capabilities and maximal 
performance characteristics of the cell. To change these 
capabilities the cell would have to acquire new biochemical 
5 reactions through some evolutionary mechanism. In so doing the 
cell would of necessity increase the range of feasible routes by 
which it could meet certain cellular demands from a set of 
environmental supplies . 

2 . 2 The Utility of Mathematics to Analyze Biochemical Reaction 

10 Networks 

The capabilities of cellular biochemical reaction networks (i) 
to produce necessary building blocks and energy requirements, and 
m (ii) to evolve the reaction pathways by which cellular 
yi production (s) is (are) realized, can be comprehensively examined 
using rigorous mathematics. Mathematical examination yields 

i y 

m results which are biochemically meaningful, serving to predict the 
performance of the biochemical reaction network. 

11 There exists one particular type of mathematical analysis of 
a_ cellular biochemical reaction networks called "convex analysis". 
3j need definition. Some of the principles of convex analysis were 
h* previously used by Schuster to find "elementary nodes", or 
!Z reactions within the biochemical reaction networks. See Schuster, 
jl S., T. Dandekar and D. A. Fell, 1999, Detection of elementary flux 

modes in biochemical networks: a promising tool for pathway 
25 analysis and metabolic engineering, Trends Biotechnology 17(2): 
53-60. See also Schuster, S. and C. Hilgetag, 1994, On elementary 
flux modes in biochemical reaction systems at steady state. J. 
Biological Systems 2(2): 165-182. Finally see Schuster, S., C. 
Hilgetag, J. H. Woods and D. A. Fell, 1996, Elementary modes of 
3 0 functioning in biochemical networks. 

Convex analysis was also previously used by Bruce Clarke to 
find "extreme currents", or reaction pathways through the 
biochemical reaction networks by which pathways the biochemical 
reaction networks succeed in processing environmental substances 
35 into the building blocks and energy requirements of the cell. See 
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Clarke, B. L., 1980, Stability of Complex Reaction Networks. 
Advances in Chemical Physics 43: 1-215. See also Clarke, B. L . , 
1981, Complete set of steady states for the general stoichiometric 
dynamical system. J". Chem. Phys. 75(10): 4970-4979. 
5 The mathematics associated with convex analysis may be used to 

determine the minimal set of biochemical pathways by which some 
particular capability of the biochemical reaction network is 
realized. These pathways satisfy both (i) mass balance constraints 
(associated with stoichiometry) and (ii) directional constraints 
10 placed on reactions (associated with thermodynamics) . 

These pathways are termed "extreme pathways", and can 
m beneficially be used to examine the functional capabilities of a 
yg biochemical reaction network. Importantly, from knowledge of these 
^ extreme pathways it is possible to determine all of the possible 
Hp combinations of reactions that need to be eliminated from the 
network to remove some particular capability (ies) of the network. 
£7 From the lists of reactions it is a simple step to determine the 
a enzymes and genes responsible for these reactions. 

^ Consider now that the elimination of these genes should then 

2jp render the biochemical reaction network incapable of reaching some 
Hj particular outcome, some particular demand (s) of the cell! 

Conversely, once it is understood what a biochemical reaction 
network is doing, and how it is doing it, then it may become 
possible to "re-engineer 11 the network, and the organism, to steer 
25 more of its output into desired channels (i.e., to make more of a 
desired reaction product) . 

Mathematical tools that permit recognition of pathways within 
biochemical reaction networks, and of the genes involved with the 
reactions within these pathways, have still further implications 
30 for the development of antibiotics to combat microbial infections. 
The tools permit recognition of how a deleterious process and 
pathway of the biochemical reaction network might be stopped, or at 
least disrupted. 

Alternatively, the same mathematical computational tools can 
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be used to improve the design and engineering of organisms for 
industrial application such as the production of bio-commodities. 
The tools permit recognition of how an beneficial process and 
pathway of the biochemical reaction network might be augmented or 
5 accentuated. 

2 . 3 Specific Prior Art Mathematical Analysis of Biochemical 
Reaction Networks 

Convex analysis has been previously used to study biochemical 
systems and to generate related sets of pathways called elementary 

10 modes and extreme currents. For a comprehensive review see the 
g paper by inventor of the present invention Schilling and his 
^0 colleagues, see Schilling, C.H., S. Schuster, B. 0. Palsson and R. 
I«s Heinrich, 1999, Metabolic pathway analysis: basic concepts and 
S3 scientific applications in the post-genomic era. Biotechnology 

11 Progress 15 (3) : 296-303. 

y= Elementary modes have been used to metabolically engineer 

;L bacteria for producing aromatic amino acid precursors at yields 
Si near the maximum theoretical yield. See Liao, J. C. , S-Y Hou and 
u Y _ p chao ^ 1996, Pathway Analysis, Engineering, and Physiological 
iij Considerations for Redirecting Central Metabolism. 
M 8 Liao, et al . , report research where all of the elementary 

modes for a reduced reaction network in Escherichia coli were 
calculated and studied to determine the optimal flux distributions 
through a central metabolism that redirected carbon flow to the 
25 pathways from aromatic amino acid production. Reactions that did 
not appear in the optimal pathways were considered indispensable, 
while those that did appear in the optimal pathways were candidates 
for over-expression. 

A similar analysis can be performed with the extreme pathways 
30 rather than elementary modes. For the precise difference between 
these two approaches see Schilling, C. H., D. Letscher and B. 0. 
Palsson, 2000, Theory for the systemic definition of metabolic 
pathways and their use in interpreting metabolic function from a 




- 5 - 



pathway-oriented perspective. Journal of Theoretical Biology 
203 (3) : 229-248 . 

Still other papers and publications discuss the complex 
analysis of biochemical reaction networks. See (i) Clarke, B. L. 
1980, Stability of Complex Reaction Networks, Advances in Chemical 
Physics A3: 1-215; (ii) Clarke, B. L. , 1981, Complete set of steady 
states for the general stoichiometric dynamical system, J. Chem. 
Phys. 75(10): 4970-4979; (iii) Edwards, J.S., R. Ramakrishna, C. H. 
Schilling and B. 0. Palsson, 1999, Metabolic flux balance analysis, 
(iv) Metabolic Engineering, S. Y. Lee and E. T. Papoutsakis, New 
York, Marcel Decker, Inc.: 13-58; (v) Liao, J. C. , S-Y Hou and Y-P 
Chao, 1996, Pathway Analysis, Engineering, and Physiological 
Considerations for Redirecting Central Metabolism, Biotechnology 
and Bioengineering 52: 129-140; (vi) Schilling, C. H., S. Schuster, 
B. O. Palsson and R. Heinrich, 1999, Metabolic pathway analysis: 
basic concepts and scientific applications in the post-genomic era. 
Biotechnology Progress 15(3): 296-303; (vii) Schuster, S., T. 
Dandekar and D. A. Fell, 1999, Detection of elementary flux modes 
in biochemical networks: a promising tool for pathway analysis and 
metabolic engineering, Trends Biotechnology 17(2): 53-60; (viii) 
Schuster, S. and C. Hilgetag, 1994, On elementary flux modes in 
biochemical reaction systems at steady state, J. Biological Systems 
2(2): 165-182; (ix) Schuster, S., C. Hilgetag, J. H. Woods and D. 
A. Fell, 1996, Elementary modes of functioning in biochemical 
networks; (x) Computation in Cellular and Molecular Biological 
Systems, R. Cuthbertson, M . Holcombe and R. Paton, London, World 
Scientific: 151-165; and (xi) Varma, A. and B. 0. Palsson, 1994. 
Metabolic Flux Balancing: Basic concepts, Scientific and Practical 
Use. Bio/Technology 12: 994-998. 
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SUMMARY OF THE INVENTION 




The present invention contemplates improvements to the 
existing mathematical method of convex analysis for purposes of 
analyzing the production of one or more selected metabolites of a 
5 biochemical reaction network so producing metabolites. In simple 
terms, the improved mathematical method of the present invention 
explains why a biochemical reaction network -- also known as a 
cellular metabolic network -- does what it does. The utility of so 
knowing is, of course, that precision understanding of the 
10 processes of nature best supports the effective, useful and safe 
q manipulation of these processes. For example, a biochemical 
process inducing disease may sometimes usefully be stopped; a 
Si biochemical process producing a valuable biochemical as an 
CQ metabolic output may sometimes usefully be enhanced. 
£§ The present invention contemplates still more than gaining 

M= knowledge of a biochemical reaction network: the present invention 
!L further contemplates an improved application of the mathematical 

ess? 

m process of convex analysis particularly for the purpose of 
N= identifying critical points in cellular metabolic networks engaging 
S in biochemical reactions. The improved analysis proceeds so that, 
by identification of some critical point, all members of the 
particular set of complex pathways of biochemical reactions leading 
to this point may be better understood. The better understanding 
this set of complex pathways again permits the better understanding 
25 of which reactions, and associated reaction pathways, can either be 
(i) disrupted so as to defeat that the cellular metabolic network 
should attain this critical point, or (ii) enhanced, so that the 
cellular metabolic network will produce more of some selected 
biochemical (s) at the critical point, and/or produce this 
30 biochemical (these biochemicals) faster. In simple possible terms, 
the improved mathematical method of the present invention not only 
explains how a biochemical reaction network, or cellular metabolic 
network, functions (as in the preceding paragraph) , but permits 



analysis of why and how the network might be (i) precluded or 
obstructed from working (as in the avoidance of disease) , or, 
alternatively (ii) made to work better and/or faster (as in the 
production of a useful biochemical product) . 

The present invention contemplates yet still more: an improved 
method of deriving a mathematical system of linear equations and 
linear inequalities representing -- i.e., serving as a mathematical 
model of -- a biochemical reaction network (also known as a 
metabolic network) in other that the derived system may be 
rigorously analyzed. Improvements to the existing mathematical 
process of convex analysis are again involved. In simplest 
possible terms, the present invention delivers more than just a 
"road map" to biochemical life processes (as in the second 
preceding paragraph) , or even a "road map" with all major "way 
points" prominently accurately identified (as in the immediately 
preceding paragraph) , but actually delivers a mathematical model of 
the biochemical reactions of a metabolic network which mathematical 
model quantitatively captures the production capabilities of the 
network . 

Clearly such knowledge of what is produced, by what processes, 
when, and to what amount (s) starts to subject life processes -- as 
are conducted by cellular metabolic networks engaging in 
biochemical reactions -- to the techniques of (bio) chemical 
engineering. One the (bio) chemical "plant" and its processes are 
understood, then a (bio) chemical engineer should, and is, able to 
approach modification of the "plant" and/or its (life) processes in 
order to obtain different, or better, results. 

1 • General Approach of the Present Invention 

The present invention provides a general framework and system 
for the identification of all the minimal sets of reactions that, 
when removed from a biochemical reaction network, will render the 
network unable to reach its particular production objectives. 

As precursor steps to the present invention, biochemical 



reaction networks are conventionally (i) constructed from genomic 
and biochemical data and (ii) described by a stoichiometric matrix. 
Together, the constraints on the directions in which reactions can 
proceed and the stoichiometric matrix correspond to a mathematical 
5 system of linear equations and linear inequalities, which system 
can be studied using convex analysis. 

In accordance with the present invention, sets of extreme 
pathways are defined that are used to represent all of the possible 
steady states which the network can achieve. By removing a single 
10 reaction in the network all of the pathways that utilize this 
^ reaction are also removed. To remove the ability of the network to 
*2 reach a particular objective, all of the extreme pathways that 
J3 reach this objective must be removed. 

It In one embodiment of the invention, a set of reactions and 

p> their associated genes are identified that will eliminate the 
ability of the target organism to generate essential biomolecules 
necessary for its growth. 

D Rather than removing one reaction at a time -- or many 

r? reactions at a time -- in order to determine the effects of the 
removal (s) on the performance of the reaction network, it is ideal 

j~f to be able to calculate all of the sets of reactions that will 
eliminate the capability of the network to achieve a particular 
production objective. In an additional embodiment of the 
invention, an algorithm is presented by which it is possible to 

25 generate the entire group of sets of reactions that can be used to 
eliminate functional properties and network capabilities of 
interest in a target organism. Collectively this group of reaction 
sets will be called "minimal deletion sets". 

These deletion sets are unique to the present invention, and 

3 0 quite valuable. They may be, for example, used to identify 
reactions which are critical to the performance of the network 
under a particular condition, but whose function is essentially 
redundant under other conditions. In addition to finding global 
deletion sets that are critical under all conditions it is then 
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possible to find environmental or condition-specific deletion sets. 

In summary, beginning with an advanced, but conventional, 
stoichiometric matrix and a set of general reaction constraints, 
the present invention shows how to generate (i) the set of extreme 
pathways along with (ii) the process for identifying a minimal 
deletion sets from these pathways. An (iii) algorithm is then 
shown by application of which algorithm all the minimal deletion 
sets in the network may be may be identified. 

These mathematical processes of the present invention, as 
discussed in greater detail hereinafter, are normally embodied in 
a computer software program. The program can be used to calculate 
the extreme pathways and minimal deletion sets of a biochemical 
yp reaction network of virtually any complexity. These valuable 
g minimal deletion sets can then be used, for example, in the 
\f subsequent development of drugs to combat infectious disease. They 
g may alternatively be used for the rational design and engineering 
s of organisms for the production of biomolecules of interest (often 

O referred to as metabolic engineering) . 

^ 2 • A Metho d of Analyzing the Production of One or More Selected 
U Metabolites o f a Biochemical Reaction Network Producing 

Metabolites 

Accordingly, in one of its aspects the present invention is 
embodied in an improvement to an existing method of analyzing the 
production of one or more selected metabolites of a biochemical 
25 reaction network that produces metabolites. 

The existing method has as inputs (i) reactions of the 
biochemical reaction network constructed from genomic and 
biochemical data, (ii) exchange fluxes on such of the produced 
metabolites as are of interest as inputs and outputs to the 
network, (iii) a stoichiometric matrix, developed from the 
reactions in consideration of the exchange fluxes, defining 
participation of each network metabolite in each reaction and 
exchange flux of the network, and (iv) a system of linear equations 
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and inequalities mathematically defining the network. From these 
inputs the existing method serves to identify deletion sets of 
reactions that, when removed from the network, eliminate the 
capability of the network to produce a selected metabolite. 
5 In particular, the improvement to this existing method in 

accordance with the present invention uses the linear equations and 
inequalities of the network to mathematically calculate a convex 
solution space called a "flux cone". The calculating produces 
"generating vectors" of this flux cone, which generating vectors 
10 are called "extreme pathways". Using these generating vectors 
(called extreme pathways) , the mathematical process of the present 
t: invention continues to determine sets of reactions that, when 
y3 deleted, diminish or eliminate capability of the network to produce 

^ y an output metabolite of interest. 

03 

15 These mathematically determined sets correspond to critical 

jO reactions of the network which, when stopped, affect the capability 

" s of the network to produce the output metabolite of interest. 
□ The method of the present invention may, after the determining 

rj of sets of reactions, optionally continue with selecting from the 

s — K 

|£ determined sets of reactions those sets that totally eliminate the 

m capability of the network to produce the output metabolite of 
interest. In this case the selected sets are called "deletion 
sets" because deletion of the reactions represented by the pathways 
of these sets suffices to totally eliminate the production of the 

25 output metabolite of interest by the network. 

The calculating of the generating vectors of the flux cone 
preferably ensues by specific mathematical manipulations within the 
more general mathematical process of convex analysis. More 
particularly, the specific preferred convex analysis of the present 

30 invention consists of calculating any of (i) a conical basis, (ii) 
a convex basis, (iii) a linear basis, or (iv) a combination of any 
of conical and convex and linear bases. 

Preferably in the method of the present invention at least 
some of the constructed reactions will have an associated 
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constraint upon the direction in which the reaction can proceed. 

The method of the present invention may be used, for example, 
to produce an output of interest which consists of one or more 
functional properties of interest in the analyzed biochemical 

5 production network. In this case the reaction sets show how these 
one or more functional properties of interest can be diminished or 
eliminated. If the output of interest consists of, for example, 
but one single functional property of interest in the analyzed 
biochemical production network then the reaction sets show how this 

10 functional property of interest can be diminished or eliminated. 

For example, that biochemical reaction network which is 

*5 analyzed by the mathematical method of the present invention can 

Jl represent a disease -producing, pathogenic, organism. The 

fy metabolite of interest will be one that is necessary for survival 

CO 

15 of the pathogenic organism. In this situation the method of the 
yB present invention using the reaction set can be directed to 
r targeting development of a drug that, by obstructing those 

□ reactions of the pathogenic organism that produce the metabolite 
'ff necessary for survival of the organism, serves to eliminate the 
4s pathogenic organism. 

O As another example, the biochemical reaction network that is 

analyzed can again represent a disease-producing, pathogenic, 
organism. However, this time the metabolite of interest can be the 
actual disease-producing, deleterious, substance that is produced 

25 by the pathogenic organism. In this situation the method of the 
present invention using the reaction set can be directed to 
targeting the development of a drug that, by obstructing those 
reactions of the pathogenic organism that produce the metabolite 
that induces disease, serves to eliminate the deleterious, disease- 

3 0 causing, function of the pathogenic organism. 

As an oppositely-directed example, when the reaction network 
analyzed is an organism producing both desired bio-molecules of 
value and un-desired bio-molecules of no value, and when the 
metabolite of interest produced by the organism is defined to be 



the un-desired and valueless bio-molecules, then the method of the 
present invention using the reaction set can be directed to 
metabolically re-engineering the organism to fail of those 
reactions that produce the particular metabolite that is un-desired 
and valueless. Thus production of un-desired valueless bio- 
molecules can be eliminated while continued production of desired 
valued bio-molecules is permitted. 

Similarly to this example, when the reaction network analyzed 
is an organism producing desired bio-molecules of value by each of 
two or more multiple -- metabolic routes, and when the 

metabolite of interest is defined to be the valued molecule as is 
produced by one only -- preferred -- route of the multiple routes 
by which the organism is capable of producing this molecule, then 
the method of the present invention using the reaction set can be 
directed to metabolically re-engineering the organism to fail of 
those reactions that produce the metabolite of interest via 
inefficient route (s), therein by eliminating production of 
metabolite via this route (these routes) nonetheless that the 
metabolite is of value. Production of the desired metabolite by 
one or more alternative one(s) of the multiple metabolic routes is 
left intact, and may even be accentuated. 

These combinations show the reaction set produced by the 
mathematical method of the present invention to be a valuable tool. 
In simple terms, the reaction set shows how to preclude, or to 
obstruct, or to accentuate individual biochemical pathways within 
the organism as lead to the production of particular metabolites. 
It is hard to ask for more than this: a complete quantitative, 
mathematical, model as to the biochemical reactions of the cell. 

3 . A Method of Identifying Critical Points in Cellular Metabolic 
Networks Engaging in Biochemical Reactions 

In another of its aspects the present invention is embodied in 
a method of (i) applying in an improved way the existing 
mathematical process of convex analysis so that a convex hull is 
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defined and spanned by unique generating, or edge, vectors, this 
hull being analyzed to derive a particular solution that is, 
mathematically, a particular point described by a flux vector lying 
within the interior of the convex hull. All this effort to get to 
5 this particular solution is for the purpose of (ii) identifying 
critical points in cellular metabolic networks engaging in 
biochemical reactions so that, by identification of some critical 
point, the particular set of complex pathways of biochemical 
reactions leading to this point may be better understood. All this 

10 effort to get to the understanding of all the biochemical reactions 
leading to a particular point is so that, (iii) by better 

O understanding this set of complex pathways, it may further be 

~ better understood which reactions, and associated pathways, can be 

yy 

ry selectively disrupted so as to defeat that the cellular metabolic 
Si network should attain this critical point. 

Iq The method of the present invention thus consists of using 

this convex hull --a mathematical construction --to represent the 
!L capabilities of a metabolic genotype. By this usage the unique 
03 generating, edge, vectors that define and that span the convex hull 
2t> represent systemically independent extreme pathways of the 
□ metabolic, life, processes of the metabolic genotype. 
^ Likewise in this usage every point in the hull is some non- 

negative combination of the unique generating, edge, vectors 
corresponding to the fact that every metabolic, life, process of 
25 the metabolic genotype is some combination of the extreme pathways 
of these metabolic processes. 

Next in the method of the present invention, the convex hull 
is mathematically solved, again by a specific application of the 
more general process of convex analysis, so as to derive a 
30 particular solution that represents a metabolic phenotype. This 
particular solution is, mathematically, a particular point 
described by a flux vector lying within the interior of the convex 
hull. 

The mathematical solving is repeated until a complete set of 
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particular solutions, corresponding to a set of flux vectors each 
lying within the convex hull, is derived. This set of solutions 
corresponds to all the pathways by which a particular metabolic 
phenotype is realized. 
5 Thus the derivation of all pathways by which the particular 

metabolic phenotype is realized is tantamount to recognition of all 
the biochemical reactions that, as part of any pathways, lead to 
the particular metabolic phenotype. Moreover, recognition of all 
biochemical reactions variously leading to the particular metabolic 

10 phenotype permits better understanding of what biochemical 
reactions of the metabolic genotype can be in particular disrupted 
so as to cause that the metabolic genotype should be unable to 

2 realize the particular solution. 

^ For example, this the method of the present invention can be 

11 employed on the genotype of a pathogenic, disease -causing, 
organism. In this case the method of the present invention 

J™ preferably continues with the development of drugs that, by 
s obstructing those biochemical reactions of the genotype of the 

pathogenic organism that lead to a particular, disease- inducing, 
p) solution of the genotype, serve to eliminate the deleterious, 
Jf disease-causing, phenotype of the pathogenic organism, 
u For example, this the method of the present invention can be 

employed on the genotype of an organism producing both (i) desired 
bio-molecules of value and (ii) undesired bio-molecules of no 
25 value. In this case the method of the present invention preferably 
continues with metabolically re-engineering the organism so as to 
obstruct those biochemical reactions of the genotype of the 
pathogenic organism that lead to that particular solution where the 
phenotype produces the undesired valueless bio-molecules, 

3 0 eliminating production of these undesired valueless bio-molecules 

while permitting continued production of desired valued bio- 
molecules . 
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4 . A Method of Analyzing a Metabolic Network 

In still another of its aspects the present invention is 
embodied in a method of analyzing a metabolic network. 

The method consists of first identifying all biochemical 
5 reactions occurring in the metabolic network, including any 
directions thereof. Then all exchange fluxes are specified, 
including any associated directional restraints attendant upon 
metabolites of the identified biochemical reactions. 

A stoichiometric matrix where each column in the matrix 
10 corresponds to a reaction, or flux, and where each row corresponds 
to a different metabolite involved in the metabolic network is next 
n created. This created stoichiometric matrix represents, in all its 
ji columns and rows, the collective biochemical reactions, being a 
Jjf form of chemical conversion, and the collective cellular transport 
3pji processes of the metabolic network, which cellular transport 
processes are how the metabolites enter and leave the metabolic 

yg 

^ network . 

s_ All directional constraints on the exchange fluxes are next 

rf combined with the created stoichiometric matrix to define the 
W metabolic network as a system of linear equations and linear 

^7 inequalities. 

Li 

U Finally the system of linear equations and linear inequalities 

that jointly define the metabolic network are analyzed by the 
mathematical process of convex analysis. 

25 The metabolic network defined as a system of linear equations 

and linear inequalities may particularly obeys the three equations 

S • v = 0 (Equation 1) 

where S refers to the stoichiometric matrix of the system and v is 
the flux vector; 



30 



Vi ;> 0, Vi 



(Equation 2) 
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where Vi corresponds to the flux value of the i th reaction; and 

Qfi s b ± s (Equation 3) 

where a t and /3 ± are ether zero of negative and positive infinity, 
respectively, based on the direction of exchange flux, and bi 
5 corresponds to the i th exchange flux. 

In this case the analyzing consists of solving the equations 
1-3 in convex space as a convex polyhedral cone in n-dimensional 
space emanating from the origin of the space. 

Every point on the convex polyhedral cone may in particular be 
represented by 

m k 

5 

L, (Equation 4) 

f7 In this case the analyzing consists of calculating the conical 

□ hull of the flux cone as representing the extreme pathways in the 
H= metabolic network . 

15 After this calculating (which is a part of the analyzing) , the 

method may further continue by determining from the calculated 
pathways critical biochemical reactions, or sets of biochemical 
reactions, that are required for the metabolic network to attain a 
particular objective or group of objectives -- as is (are) 

2 0 represented by one or more particular points on the flux cone. 



These and other aspects and attributes of the present 
invention will become increasingly clear upon reference to the 
following drawings and accompanying specification. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Referring particularly to the drawings for the purpose of 
illustration only and not to limit the scope of the invention in 
any way, these illustrations follow: 
5 Figure 1 is a flow diagram illustrating one procedure of the 

present invention for determining one set, or the entire 
collection, of minimal deletion sets that eliminate particular 
production capabilities of interest in a metabolic network. 

Figure 2 is a geometric representation of the flux cone of the 
10 mathematical method of the present invention shown in three - 
dimensions where the entire unbounded flux cone is spanned by the 
generating vectors representing the capabilities of a metabolic 
genotype . 

m Figure 3, consisting of Figures 3a-3c, are respectively a 

tj? graph of an exemplary metabolic reactions scheme, a legend to the 
2 mathematical representation of the scheme, and table of the extreme 
s pathways as collectively define a metabolic genotype and phenotype 

S in the context of convex analysis, where, more particularly, 
M= Figure 3 (a) shows a hypothetical reaction network comprised 

25 of a series of internal and exchange fluxes, functioning to 
|I generate appropriate ratios of metabolites C, D, and E for 

incorporation into biomass represented as the GRO metabolite; 

Figure 3 (b) shows a mathematical translation of the reaction 

network into the steady- state mass balances and constraints placed 
25 on all fluxes that define the solution domain; and 

Figure 3(c) is a list of the ten (10) extreme pathways 

calculated for the network shown in vector format indicating the 

relative flux activity of each of the fluxes in the pathways. 

Figure 4 is a table 1 showing the set of extreme pathways 

30 *pl, ,pl0) for the reaction scheme shown in Figure 3. 

Figure 5 is a graph showing normalized values for the 

objective flux (b z ) for all single and double deletion combinations 

of the network described in Figure 3 . 
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DESCRIPTION OF THE PREFERRED EMBODIMENT 

Although specific embodiments of the invention will now be 
described with reference to the drawings, it should be understood 
that such embodiments are by way of example only and are merely 
5 illustrative of but a small number of the many possible specific 
embodiments to which the principles of the invention may be 
applied. Various changes and modifications obvious to one skilled 
in the art to which the invention pertains are deemed to be within 
the spirit, scope and contemplation of the invention as further 
10 defined in the appended claims. 

The present invention relates to systems and methods for (i) 
%D identifying deletion sets in a target organism's metabolic network 
r; (ii) based on a pathway analysis (iii) mathematically derived from 
03 (iv) a list of reactions associated with convexity, which reactions 
fS are a partial or a complete representation of reactions in a cell. 
14 A stoichiometric matrix is used to describe the connectivity 

1^ amongst all the metabolites in a metabolic network of any 
S complexity including that for an entire organism. Additional 
constraints on the direction in which reactions and fluxes can 
23 proceed in the network serve to complete the mathematical 
M= description of the reaction network. From this description 
existing algorithms can be used to calculate the set of extreme 
pathways that form the conical hull of the solution space, 
corresponding to all feasible steady state flux distributions. The 
25 details of performing the construction of a stoichiometric matrix 
and the formulation of the constraints imposed on a biochemical 
reaction network have been previously disclosed along with an 
algorithm for calculating the set of extreme pathways (Schilling, 
2 000) . From these pathways it is possible to identify sets of 
3 0 reactions that upon removal from the network eliminate the ability 
of the network to achieve selected objective of interest. 
Moreover, it is possible to calculate all of these deletion sets 
using one algorithm that operates on the set of extreme pathways. 
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These minimal deletion sets may be useful for the identification 
and development of potential protein and genetic targets for anti- 
microbial drugs and the engineering of microbial strains for 
bioproduction purposes. 

5 It should be noted that the systems and methods described 

herein can be implemented on any conventional host computer system, 
such as those based on Intel® microprocessors and running Microsoft 
Windows operating systems. Other systems, such as those using the 
UNIX or LINUX operating system and based on IBM®, DEC® or Motorola® 

10 microprocessors are also contemplated. The systems and methods 
described herein can also be implemented to run on client -server 

^ systems and wide-area networks, such as the Internet. 

yy Software to implement the system can be written in any well- 

?! known computer language, such as Java, C, C++, Visual Basic, 

iii FORTRAN or COBOL and compiled using any well-known compatible 

HJ" compiler. 

The software of the invention normally runs from instructions 
a stored in a memory on the host computer system. Such a memory can 

Ji be a hard disk, Random Access Memory, Read Only Memory and Flash 
W Memory. Other types of memories are also contemplated to function 
^ within the scope of the invention. 

y[ The process 10 for determining all minimal deletion sets for 

any metabolic network of interest is shown in Figure 1. Beginning 
at a start state 12, the process 10 then moves to a state 14 to 

25 gather all of the known biochemical information pertaining to the 
reactions in the network. This generates a list of all the 
chemical reactions occurring in the network. This provides the 
stoichiometry of each chemical reaction proposed to occur in the 
network along with information regarding the irreversible and 

30 reversible nature of the reactions. The stoichiometry of each 
reaction provides the molecular ratios in which reactants are 
converted into products. These reactions include all transport 
reactions, enzymatic reactions, and diffusion processes. For an 
entire organism these reactions are catalyzed by gene products 
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whose potential presence in the organism can be inferred from the 
annotated genome sequence of an organism along with additional 
biochemical information. Thus the reactions determined in state 14 
represent all of the physical -chemical conversions that can 
5 reasonably occur in the system. 

A theoretical system boundary can be drawn around all of the 
physically occurring reactions and associated substrates and 
products defined in state 14. These reactions contained within the 
system boundary will be collectively referred to as the internal 
10 fluxes of the system. The entire group of substrates and products 
are collectively referred to as the metabolites of the system. The 
p process 10 then moves to a state 16 wherein all of the exchange 
fluxes are specified in the system. These exchange fluxes will 
JT| constitute the presence of potential sources and/or sinks on 
tB individual metabolites. Therefore if a particular metabolite is to 
r Z be allowed to enter the system or exit the system an exchange flux 
M» is created to allow for the passage of the metabolite across the 
theoretical system boundary. Conceptually these fluxes can be 
m thought of as the input and outputs of the system and can be 
26 defined by the researcher to simulate environmental conditions of 
interest. 

The information obtained in state 14 and 16 regarding the 
reactions and metabolites participating in the system must be 
translate into a mathematical format to facilitate the remainder of 

25 the process 10. While reversible reactions may be represented as 
one reaction, for ease of calculation all reactions that are 
reversible will be decomposed into a forward reaction and a 
backward reaction that are constrained to only proceed in one 
positive direction. Therefore all internal fluxes will be 

3 0 constrained to take on only positive values representing their 
activity levels. Exchange fluxes that serve as both input and 
outputs for a particular metabolite may also be decomposed in a 
similar fashion or remain bi-directional, thus capable of taking on 
any real value (positive and negative) . All of the information 
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regarding the internal fluxes, exchange fluxes, and their 
substrates and products can be represented in a matrix format 
typically referred to as a stoichiometric matrix. Each column in 
the matrix corresponds to a given reaction or flux, and each row 
5 corresponds to the different metabolites involved in the system. 
Thus, a given position in the matrix describes the stoichiometric 
participation of a metabolite (listed in the given row) in a 
particular flux of interest (listed in the given column) . Together 
all of the columns of the stoichiometric matrix represent all of 

10 the chemical conversions and cellular transport processes that are 
determined to be present in the network. This includes all 

□ internal fluxes operating within the system and exchange fluxes 
operating on the system as inputs and outputs. Thus, the process 

fil 10 moves to a state 18 in order to formulate all of the fluxes into 

Jll a stoichiometric matrix. 

in Next, the process 10 moves to a state 20 wherein information 

regarding the directional constraints placed on the fluxes is 

m combined with the stoichiometric matrix to fully define the 

£0 metabolic network as a system of linear equations and linear 
T6 inequalities to be analyzed using principles of convex analysis. 

Q In studying metabolic networks the principle of conservation of 

^ mass is applied in the generation of transient mass balances around 
each metabolite in the system. Each mass balance constitutes a 
differential equation describing the change in concentration of the 

25 metabolite as function of the activities of the fluxes that serve 
to generate and dissipate the metabolite. Structural aspects of a 
metabolic network are time invariant allowing the functioning of 
the network to be placed into a steady state. Eliminating all the 
time derivatives obtained from dynamic mass balances around every 

30 metabolite in the metabolic system, yields the system of linear 
equations represented in matrix notation, 



S • v = 0 



(Equation 1) 
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where S refers to the stoichiometric matrix of the system, v is the 
flux vector. This equation simply states that over long times, the 
formation fluxes of a metabolite must be balanced by the 
degradation fluxes. Otherwise, significant amounts of the 
5 metabolite will accumulate inside the metabolic network. Applying 
equation 1 to our system we let S now represent the stoichiometric 
matrix constructed in state 18. 

To completely describe the metabolic system it is necessary to 
include the constraints on the possible directions of the internal 
10 and exchange fluxes. The constraints on the internal fluxes is 
rather straightforward as all fluxes must be non-negative yielding: 

JJ Vi ;> 0, Vi (Equation 2) 

fQ where ^ corresponds to the flux value of the i th reaction. The 
mf Z constraints on the exchange fluxes depend on the status of the 
5 determined source or sink on the associated metabolite, or 
3 similarly on the input and output status of the metabolite. This 

can be expressed in Equation 3 below where a ± and are either zero 
M or negative and positive infinity, respectively, based on the 
2 direction of the exchange flux and b ± corresponds to the i th exchange 
P flux. 

«i < bi 5 (Equation 3) 

Under the existence of a source (input ) only a t is set to negative 
infinity and ft ± is set to zero, whereas if only a sink (output) 
exists on the metabolite is set to zero and f$ ± is set to positive 

25 infinity. If both a source and sink are present then the exchange 
flux is bi-directional with ot ± set to negative infinity and is 
set to positive infinity leaving the exchange flux unconstrained. 

Together equations 1-3 describe the metabolic system under 
steady-state conditions as a system of linear equalities and linear 

30 inequalities. The presence of linear inequalities limits the use 
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of traditional concepts of linear algebra, and necessitates the use 
of convex analysis, which is capable of treating systems of linear 
inequalities. The set of solutions to any system of linear 
inequalities included that described above is a convex space. This 
5 convex solution corresponds geometrically to a convex polyhedral 
cone in n-dimensional space {R n ) emanating from the origin for all 
metabolic systems modeled as described herein. We refer to this 
convex cone generally as the flux space and more specifically as 
the steady-state flux cone (C) . Within this flux cone lies all of 
10 the possible solutions and hence flux distributions under which the 
system can operate. Since every solution or operating mode of the 

□ system is contained within the flux space, it logically follows 
^ that the entire flux space represents the capabilities of the given 
fjj metabolic network. Thus the flux space clearly defines what a 
§i network can and cannot do . 

^ Convex cones are described by the extremal rays (or generating 

vectors) that correspond to the edges of the cone, (being half- 
;L lines emanating form the origin) . These extremal rays are said to 
Eg generate the cone (forming the conical hull) and cannot be 
ft) decomposed into a non- trivial convex combination of any other 

□ vectors residing in the flux cone. Here in the context of 
H= metabolic systems the term extreme pathways is used to denote the 

extreme rays of a polyhedral cone as each ray corresponds to a 
particular pathway or active set of fluxes which satisfies the 

25 steady-state mass balance constraints and inequalities of equations 
1-3. Extreme pathways are denoted by the vector p t and the total 
number of extreme pathways needed to generate the flux cone for a 
system will be denoted by k. Every point within the cone (C) can 
be written as a non-negative convex linear combination of the 

3 0 extreme pathways as shown below: 

k 

C = v: v = Y, w i p i ,co i >0 V i 



(Equation 4) 



Thus the set of extreme pathways is analogous to a basis/coordinate 
system that can be used to describe a position in space. These 
pathways are said to conically span or generate the set of all 
pathways as any pathway and/or distribution of fluxes can be 
written as a non-negative linear combination of the p/s. The 
pathway vector w corresponds to the coordinate vector relative to 
the set of extreme pathways. It provides the weight given to each 
pathway in a particular flux distribution (v) . 

The calculation of the unique set of extreme pathways for a 
metabolic network, described by equations 1-3, moves the process 10 
to a state 22. These pathways can be calculated using any 
algorithm capable of generating the conical hull of a convex 
polyhedral cone. One such algorithm has been previously disclosed 
(Schilling, 2000) . 

The set of extreme pathways geometrically represent the edges 
of the flux cone as depicted in Figure 2, where the flux cone is 
the region of all admissible flux distributions that a metabolic 
network can display under the assumption of steady state. These 
pathways define the structural capabilities of the network and can 
be used to interpret every possible flux distribution or metabolic 
phenotype/functional modality of a metabolic network. Thus they 
define what a network can and cannot do; what building blocks can 
be produced; how efficient can the network generate biomass or 
extract energy from substrates; what are the redundancies in the 
network . 

Of particular interest is the ability to use these pathways to 
identify critical reactions or sets of reactions that are required 
for the network to reach a particular objective (s) . These lethal 
or minimal sets of reactions and the genes coding for the gene 
products of these reactions are herein referred to as minimal 
deletion sets. Furthermore, these pathways can be used to 
calculate every possible deletion set in a network for a given set 
of input and output conditions. Any combination of genetic 
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deletions that would be lethal to an organism must contain a subset 
of genes that form a minimal deletion set. Therefore if the 
combined loss of gene A and B is lethal to an organism, where A and 
B alone are not essential genes, then the set is a minimal deletion 

5 set, making every superset of this also a deletion set (such as the 
loss of gene A, B and C) . 

When considering the ability of a metabolic network to produce 
a defined set of metabolic precursors used to generate all of the 
components of the biomass, these deletion sets correspond to 

10 targets or combined targets for antimicrobial therapeutics. When 
considering the specific production of a metabolite (i.e. amino 

n acid) these sets correspond to reactions that can be deleted to 

^ direct the flow of metabolic resources in the cell for the 

51 metabolic engineering of an organism. 

; is? 

J§ The identification of the complete group of minimal deletion 

r ^ sets for a metabolic network moves process 10 to a state 24 . The 
M= procedure for determining a minimal deletion set follows a simple 
rationale. When an internal flux is eliminated from the network or 

f 5 

m identically if the flux value is forced to zero it is seen from 
2t) Equation 4 that each of the pathways that utilize the particular 
U flux are forced to be weighted to zero, otherwise non-zero values 
would exist in v corresponding to the eliminated reaction. 
Therefore all pathways that use this flux can be eliminated for 
further analysis of the situation. After removing the appropriate 

2 5 pathways, if there are no pathways capable of producing a desired 

output or objective then there exists no combination of pathways 
that can generate the particular objective. To calculate the 
complete group of minimal deletion sets an algorithm can be 
designed to check for all of the possible combinations of deletions 

3 0 that will eliminate all of the pathways capable of producing the 

desired objective. 

(XXX- yes, insert details of the computational algorithm) 
Together the formulation of the stoichiometric matrix and 

reaction constraints for a metabolic network, followed by the 
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calculation of the set of extreme pathways and lastly the 
identification of the entire collection of minimal deletion sets 
terminates the process 10 at an end state 26. 

Thus by mathematically describing the reactions deemed present 
5 within a cellular metabolic network and adding or removing 
constraints on internal and exchange fluxes (inputs) in the network 
it is possible to (1) simulate a genetic deletion event or complete 
enzymatic inhibition and (2) simulate the presence or absence of 
metabolic resources available to the network perhaps to mimic its 

10 in vivo environment. Based on the resulting mathematical 
description of the network the unique set of extreme pathways, 

q which define the complete range of metabolic phenotypes and 
*S production capabilities of the network, can be calculated. From 

i IT'S 

11 these pathways particular fluxes and/or sets of fluxes can be 

i y 

£S deleted from the network to assess the ability of the network to 
\T produce the constituents of the biomass or any particular metabolic 
objective. If the removal of a set of internal fluxes from the 
^ network eliminates the ability of the network to produce the 
m desired metabolic objective such as precursors for growth, and the 
2€ cell can not obtain these precursors from its environment, then the 
removal of this set of fluxes constitutes a minimal deletion set 
that has potential as an antimicrobial drug target or combination 
of targets. All of these minimal deletion sets can be calculated 
based on the process described herein and the appropriate genetic 
25 targets identified from the network composition. 

In addition to using these deletion sets for the 
identification of antimicrobial drug targets, these deletion sets 
can be used to design and engineer cells to have desired metabolic 
characteristics and capabilities that can be assessed using a 
30 pathway analysis. This has clear potential in the metabolic 
engineering of bacterial strains to be competitively co-cultured. 
Cells can thus be engineered to require specific growth 
environments and substrate dependencies. Opposing strains can even 
be designed to be co-dependent through the elimination of genes as 
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- 27 - 

guided by the process described herein of identifying minimal 
deletion sets. 

It is also possible to calculate all deletion sets that will 
render the network to operate sub-optimally below a certain 
5 threshold from the extreme pathways, not just deletion sets that 
result in complete functional losses. Following the same procedure 
disclosed herein with only a slight modification does this 
calculation. Instead of looking for deletion sets that eliminate 
all of the pathways that generate the particular objective, only 
10 deletion set that eliminate pathways that generate the particular 

objective above a specified yield are calculated. 
^ As a stoichiometric matrix or a connectivity matrix can be 

ji generated to describe biological networks in other areas such as 
tfS signal transduction, the identification of minimal deletion sets is 

if" 3 

iM applicable for identify drug targets beyond metabolic networks, 
extending even into the identification of critical links or sets of 
rf points in non-biological networks. 

Jr? 2 . Example 1: Minimal deletion sets in a hypothetical reaction 

: i s 

y* network 

^6 Using the process disclosed in Figure 1, the entire collection 

ri of minimal deletion sets has been calculated for a hypothetical 
reaction network to completely illustrate the approach. The 
hypothetical reaction scheme considered is shown in Figure 3 (A) . 
This network is comprised of 6 metabolites, 8 internal fluxes, and 

25 5 exchange fluxes. Note that internal fluxes v 3 and v 4 may 
represent two separate reactions or be the result of a 
decomposition of a reversible reaction into both a forward and a 
reverse reaction. We will only allow metabolites A and C to enter 
the system as available substrates, and allow metabolite GRO to 

30 exit. GRO is a metabolite used to represent the result of a 
systemic demand (i.e. a growth flux, v z ) of one mole of C and D and 
2 moles of E (analogous to biosynthetic precursors used to generate 
biomass or actual components of the biomass) . 
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After determining the reactions in the network and the 
exchange fluxes, all of the governing mass balance equations are 
assembled. These equations are provided in Figure 3(B) and can be 
used to generate the stoichiometric matrix for the system. In 
5 addition all of the directional constraints on the reactions, due 
to systems specific conditions on the input/output of various 
substrates and products and the thermodynamics of the internal 
reactions, are also provided in Figure 3(B). The next stage in 
determining the minimal deletion sets for a metabolic network is 
10 the identification of the set of extreme pathways for the network. 

For the system provided in Figure 3 (A) the calculated extreme 
p pathways are provided in Figure 3(C) . Pathway 1 and 2 use both A 
^ and C to reach the growth objective while pathway 3 and 4 only use 
/ HJ C/ and 5 ' 6/ 7 ' and 8 use only metabolite A as their input. 
tB Pathway 9 and 10 correspond to internal cycles in the network and 
J are ignored from further consideration, as they do not add 
U functionality to the network. Each pathway satisfies the mass 
JL^ balances and flux constraints placed on the network in Figure 3(B) . 
03 As previously mentioned when a reaction is eliminated from the 

M network or identically if the flux value is forced to zero we can 
p see from Equation 3 that each of the pathways that utilize the 
M= reaction are forced to be weighted to zero, otherwise non-zero 
values would exist in v corresponding to the eliminated reaction. 
Therefore all pathways that use this flux can be eliminated for 
25 further analysis of the situation. As an example consider the 
elimination of flux v x from the network in Figure 3 (A) . We can 
seen from Table 1 that this would remove pathways 1, 2 and 5 
through 8, leaving only pathway 3 and 4 as feasible pathways that 
can operate in the network. This is logical as the elimination of 
30 this particular flux would not allow metabolite A to be utilized, 
and subsequently all the pathways that are eliminated represent 
those that utilize metabolite A. 

Now if we consider the case that only metabolite C is 
available to the network by forcing Jb 2 to equal zero, we are left 
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once again with only pathway 3 and 4. The minimal deletion sets 
for growth on any particular substrate or combination of substrates 
corresponds to all of the sets of reactions that when deleted 
eliminate all of the pathways that are available to produce biomass 
5 (in our case demonstrated by a positive flux level on flux b 3 ) . 
Thus for the situation of growth on C we can see that the 
elimination of fluxes v 2 through v 4 does not affect the remaining 
pathways (p 3 and p 4 ) . Eliminating reaction v 5 and v 8 individually 
eliminates both of the remaining pathways while eliminating v 6 and 

10 v 7 together eliminate both of the pathways. Thus the minimal 
deletion sets for growth on substrate C is simply { (v 5 ) , (v 8 ) , {v 6 

q /VV)}- We can perform similar calculations for growth on substrate 
yCI A alone by focusing on pathways 5 through 8 to reveal that the 
minimal deletion sets under this condition are { (v 2 ) , (v 2 ) , (v a ) , (v 3 

11 r v s) * (^s,^?)}- For growth in the presence of both substrates we 
would have to consider all eight pathways. In this case the 

l~ minimal deletion sets are { {v 8 ) , (v 2 ,v 5 ), (v 3/ v s ), (v ff/ v 7 )}. Thus we 
examine a wide range of substrate availability conditions using 
this approach and determine the minimal deletion sets for each 
particular condition. This process of determining the minimal 

Jl deletion sets from a set of vectors can be easily automated on a 

yi computer to handle much larger reaction networks and sets of 
pathways representative of realistic cellular networks. 

Linear optimization techniques such as those used in flux 

25 balance analysis (Varma, 1994, Edwards, 1999) can also be used to 
determine the results of all the single and double deletion 
combinations to confirm our results. To consider the case that 
both A and C are available a linear programming problem can be set 
up to optimize for b 3 while constraining the exchange fluxes for A 

3 0 and C not to exceed the input of more than one mole of A and one 
mole of C into the system. The maximum value for b 3 in this case 
would equal 0.5. From this we can compare the maximum value of b 3 
for each single deletion and double deletion case to the u wild 
type" case. Figure 4 displays the results of this type of 



- 30 - 



analysis. There is no single reaction that is critical, but the 
combinations of (v 1/ v 5 ) f (v 3 ,v 5 ) , and {v 6 ,v 7 ) are lethal. These 
results are in exact agreement with minimal deletion sets 
calculated from the extreme pathways . 

In this example the demand considered on the network was that 
of the GRO flux, but we could also consider a simpler case in which 
the demand is only one metabolite such as D or E representative of 
a particular cellular demand. The process is identical with the 
exception of the inequality constraints that are placed on the 
exchange fluxes before the extreme pathways are calculated. 

In accordance with the preceding explanation, variations and 
adaptations of the mathematically based method of analyzing 
biochemical reaction network, also know as metabolic networks, in 
accordance with the present invention will suggest themselves to a 
practitioner of the mathematical arts. 

In accordance with these and other possible variations and 
adaptations of the present invention, the scope of the invention 
should be determined in accordance with the following claims, only, 
and not solely in accordance with that embodiment within which the 
invention has been taught. 



